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FINANCIAL INSTRUMENT PORTFOLIO CREDIT EXPOSURE EVALUATION 

The present invention relates to a method for use in the evaluation 
of the credit exposure of a portfolio of financial instruments. 

Financial instruments such as forwards, options and other derivatives 
can be highly complex, making an assessment of their future value 
over time a complicated and highly specialised area. 

It will be appreciated that the value of such instruments at any one 
future point in time will depend not only on the "internal" nature of 
the instrument itself, but also on "external" market factors, such as 
for example future exchange rates, asset prices and interest rates. 

The "credit exposure" of a portfolio of financial instruments is 
defined - for a future point in time - as the statistical upper- 
bound (usually 97.5% or 99%) of the possible values which might be 
attained by that portfolio at that future time. The construction of 
future values may also incorporate the effect of legal contracts 
which reduce credit exposure in the event of client default such as 
collateral agreements and close-out netting agreements. 

A reliable estimate of credit exposure (statistical upper-bound of 
future value) for a portfolio of such financial instruments can be 
invaluable to a financial institution. In particular, a profile of 
credit exposure at many future points in time is often used when 
assessing whether to enter into a financial contract with a client. 
For example, a "Tier 1" bank, such as ABN-AMRO may be considering 
whether to trade a particular financial instrument with a client, say 
Shell. ABN-AMRO need to be sure however that the addition of that 
financial instrument to Shell's existing portfolio does not create 
inappropriate future credit exposures: (a) that Shell will be able to 
meet its worst case obligations on the portfolio and (b) that the 
level of credit exposure is proportional to the risk of Shell being 
put into receivership. 
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In this connection; mathematical systems are currently available for 
analysing a portfolio of financial instruments with a view to 
predicting, at a future point in time, the statistical upper-bound of 
the value of the portfolio. The statistical upper bound of the 
portfolio can be seen as the worst case scenario for what could be 
lost, i.e. the maximum that may be lost if the counterparty cannot 
cover that maximum value due to default. 

One such method for assessing the statistical upper-bound utilises 
the Monte Carlo technique. The method requires that each financial 
instrument can be "valued" using a mathematical formula which gives 
the future value of the instrument from hypothetical future values of 
underlying market variables (asset prices, exchange rates, interest 
rates) . The value of a portfolio at some arbitrary future date can 
hence be expressed as the combination of valuation functions of the 
constituent financial instruments. This results in a very large 
mathematical function, say P() , which depends on a very large number 
of unknown market variables. 

The Monte Carlo method establishes the upper-bound of portfolio value 
by generating future random values for all the variables that the 
portfolio depends upon and applying the portfolio valuation function, 
P(), to those variables. If this is done often enough (10,000 is 
typical) then one can estimate the upper-bound at a given point in 
time by ordering the outcomes an picking the valuation that occurs at 
the appropriate rank. 

In other words, at each designated point in time, the method will 
carry out thousands of calculations to evaluate all potential values 
of the portfolio at that point. The reliability of the estimate will 
increase for an increased number of calculations, but this can be 
prohibitive since the speed of securing a result is directly related 
to the number of calculations made. The Monte Carlo method is 
laborious and hence relatively slow. 

Since the portfolio valuation function can be so complex, there is no 
guaranteed shortcut for carrying out all this multi-dimensional 
statistics without simulation. 
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An object of the present invention is therefore to seek to provide an 
evaluation method that alleviates such problems of existing methods. 

According to a first aspect of the present invention there is 
provided a method for evaluating the credit exposure of a portfolio 
of one or more financial instruments, the method comprising 

establishing a deal object for the or each financial 
instrument, the deal object comprising a representation of said 
financial instrument and a valuation function for representing how 
the value of the financial instrument is related to underlying market 
variables . 

establishing one or more risk factor models, the or each risk 
factor model representing an underlying financial market variable 
which may affect the value of one or more of said financial 
instruments; 

establishing a deal parabolic function for representing the or 
each deal object valuation function, by operation of said the or each 
deal object valuation function on said the or each risk factor model 
to which it is sensitive, each deal parabolic function giving its 
value at a particular instant under a particular state of market risk 
factors; 

summing the coefficients of each said deal parabolic function 
established at a same instant from the one or more deal objects 
represented in the portfolio in order to build a portfolio parabolic 
function which approximates the overall portfolio value for that 
instant; 

wherein establishing each said deal parabolic function involves 
evaluating a plurality of coordinates calculated from said deal 
object valuation function and then applying a parabolic curve, 
surface or multi- dimensional surface (manifold) to fit said 
coordinates, the parabolic curve or surface then representing that 
deal object valuation function. 

Preferably, three or more coordinates are calculated and plotted to 
which the parabolic curve or quadratic polynomial is to be applied. 
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The evaluation method of the present invention is not hence based 
around using random elements, but instead simplifies the problem of 
valuing complex functions so that an answer can be determined by 
linear algebra and established properties of quadratic forms. 

In this regard, the evaluation method of the present invention 
involves the creation of a simplified approximation of each deal 
object valuation function, a hyper-dimensional surface or manifold. 
This manifold is preferably quadratic (in any one dimension only 
depends on three parameters) . 

In preferred embodiments, the or each deal object valuation function 
is established using a deal object, each deal object operating to 
select an appropriate skeleton or template representation from a 
store of such representations of financial instruments, and to 
populate the skeleton or template with data which will define the 
nature of the financial instrument. 

Conveniently, the or each deal object further operates to identify to 
which risk factor models the valuation function of that deal object 
is sensitive. 

Preferably, the or each deal object further operates to value the 
deal represented therein in relation to the risk factor models that 
have been identified as being appropriate to the valuation function 
of that deal object. For a single future point in time several 
valuations are performed corresponding to different state of the 
contingent risk factor models . 

In preferred embodiments, the deal object operates to apply an 
"optimal fit" parabolic curve to the coordinates evaluated from the 
deal valuation function, the optimal fit curve being chosen on the 
basis that it most readily passes through the coordinates plotted 
with minimum deflection between coordinates. The exact position and 
number of the valuations conducted depends on the type of deal 
object. Where three points per risk factor are evaluated, an exact 
parabola fit can be established. Where more than three points per 
risk factor are evaluated the technique of quadratic least squares is 
employed to derive the parabola coefficients. 
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Preferably, the coordinates to which the parabola is to be applied 
are chosen from each of the limits of the deal object valuation 
function (defined by the extremes of the relevant risk factors) and 
at its central region. 

Deal objects can be of various types depending on the type of 
financial instrument represented. For example a Foreign Exchange 
Option will be represented by an FXOptionDeal object. This object 
will have parameters which store all the details (terms & conditions) 
of that particular deal e.g whether it is a 'call' or 'put' , whether 
it is bought or sold, the option strike amount, the underlying 
amount, the trade date, and the exercise date. 

If a portfolio contains 1000 foreign exchange options then it will be 
represented by 1000 deal objects, all with different data fields, but 
with the same type. The structure of the data they hold and the code 
that is operable on them will be the same. But of course, portfolios 
are likely to contain many types of financial instrument. In the 
current system there are around 50 different types of deal objects. 
Templates or skeleton representations of such deal objects are stored 
in a suitable store. 

Preferably, the or each risk factor model is established using a risk 
factor object, each risk factor object operating to select an 
appropriate skeleton or template representation from a store of such 
representations of risk factors and to populate the skeleton or 
template with data which will define the nature of the risk factor 
model . 

The value of a deal at any point in time depends on certain market 
variables which are termed risk factors. For example the value of a 
GBP : USDOLLAR option deal is in essence determined by the behavior of 
(a) the GBP: USD exchange rate and (b) the USD interest rate (if we 
are valuing the option in USD) . The option will only be affected by 
other market variables in as much as they affect these two 
"under lyings" . The same is true of the portfolio as a whole . A 
counterparty portfolio might contain 1000 deals and its value 
dependent on 100 risk factors. Templates or skeleton representations 
for such risk factors are stored* in a Risk Factor List. 
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Conveniently, the method further comprises compiling a database or 
store of financial and market data. 

Preferably, the method comprises the use of techniques of linear 
algebra and quadratic forms to evaluate the statistical upper-bound 
of the parabolic function for the overall portfolio, and hence, an 
approximation for the credit exposure of the portfolio. 

In preferred embodiments, said techniques comprise one or more of 
aggregation , decomposition of the covariance matrix , transformation 
onto orthogonal variables, the application of PCA (Principal 
Component Analysis) on weighted factors, the transforming away of 
cross terms, and the calculation of distribution moments and 
percentiles for the approximation. 

According to a further aspect of the present invention there is 
provided a system for evaluating the exposure of a portfolio of one 
or more financial instruments, the system comprising: - 

means for establishing a deal object for the or each financial 
instrument, the deal object comprising a representation of said 
financial instrument and a valuation function for representing how 
the value of the financial instrument is related to underlying market 
variables; 

means for establishing one or more risk factor models, the or 
each risk factor model representing an underlying financial market 
risk factor which may affect the value of one or more of said deal 
objects; 

means for establishing a deal parabolic function for each deal 
object, comprising means for operation of the deal object valuation 
function on said the or each risk factor model to which it is 
sensitive; and 

means for summing each said deal parabolic function established 
at a same instant from the deal objects represented in the portfolio 
in order to build an overall portfolio value for that instant; 

wherein said means for establishing each said deal parabolic 
function comprises means for making an approximation of said each 
deal object valuation function, involving the evaluation of a 
plurality of coordinates calculated from said deal object valuation 
function and then applying a parabolic curve, surface, hyper surface 
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(manifold) to fit said coordinates, the parabolic curve or quadratic 
polynomial than representing that deal valuation function. 

A system as defined above, wherein the system is a computer system. 

The present invention encompasses a computer program comprising 
program instructions for carrying out the above method a computer 
program as defined above, the program being embodied on a record 
medium, computer memory, read-only memory, or electrical carrier 
signal . 

An example of the present invention will now be described with 
reference to the attached drawing, Figure 1, which shows 
diagrammatically the method according to the present invention. 

It will be appreciated that the method that forms the present 
invention can be implemented on a computer based machine. As such the 
elements represented in Figure 1 correspond to computer system (a 
computer-based machine including computer hardware and software) . 

In this regard a computer system performing the method of the present 
invention, is used to calculate the upper-bound of the possible 
future values of a portfolio of financial instruments . The 
calculation can be equated to a maximum exposure of that portfolio 
and hence can be used by a financial institution to evaluate whether 
a particular portfolio is appropriate for a particular counterparty. 

As a first step in the method, predetermined mathematical 
representations of the financial instruments to form the portfolio 
are called up from an electronic memory store and populated with data 
from a database to thereby define the nature of each of the 
respective financial instruments. 

Hence, for each type of financial instrument to the portfolio there 
exists a electronic basic skeleton or template representation. Such a 
skeleton is matched to the financial instrument in question, and then 
populated with data from a data base, to define the financial 
instrument more precisely. As an example, there is a skeleton 
representation for certain foreign exchange related financial 
instruments. The representation will for example have parameters 
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which store details of whether the instrument is a "call" or "put", 
whether it is bought or sold, a strike amount, an underlying amount, 
a trade date and an exercise date. Such data is populated into the 
skeleton function to define the nature of the instrument. 

The defined representations of the financial instruments are termed 
"deal objects". The deal objects have three main functions: 

The first function of each deal object is to value itself using its 
deal valuation function for a given scenario of market variables. 
The valuation function is an established mathematical function which 
is different for each deal type. The value of the function depends 
on the values of market variables represented on the risk factor 
objects and also the representation of the financial instrument 
contained in the deal object. 

The second function of each deal object is to establish to which 
external or market factors the valuation function of that deal object 
is sensitive. For example, for a foreign exchange financial 
instrument, factors such as the USDollar/Sterling rate can have 
significant effects on their valuation. As such, the second function 
of each deal object is to identify which external factors (termed 
risk factors) are appropriate to the deal valuation function. In 
this connection, the future value of a financial instrument (or 
portfolio of instruments) will be affected by one or more inherently 
random market variables - i.e. the risk factors. Thus when making a 
future valuation of a particular instrument, one needs to take into 
account the effects of such underlying factors. 

To this end, risk factor objects operate within the method to 
evaluate the effects of such external risk factors. The risk factor 
objects operate in a similar way to the deal objects, namely for each 
type of risk factor, such as asset prices or interest rates, there is 
an electronic skeleton or template in the form of a mathematical 
representation. These representations are populated with data from a 
database to define more precisely what that risk factor is and 
thereby form one or more risk factor models. The risk factor objects 
are hence a statistical description of the market variable in 
question, say UK interest rates, together with a particular state of 
that variable. The "state" of each risk factor is used throughout 
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the method in order to represent different future market scenarios to 
the deal valuation functions. 

The third function of each deal object is to approximate its 
valuation function with a parabolic function. The parabolic function 
for each deal object is determined using the deal valuation function 
and the relevant risk model or models, and is represented by the 
formula : - 

U(x(t)> 0 = <*(0 + ^ bi(t). Xi (t) + 1 a/OjcffljCjft) 

Where : - 

UfefOtO is tne parabolic approximation of the valuation 

function. 

x(t) is multivariate normal vector of risk factors. 

a, b f C are time -dependent constants . 

In this regard, each deal object is able to select how its deal 
valuation function fits a parabola in the "optimal" way such that its 
behavior at the extremes and any internal nuances are captured. 
Fitting the "optimal fit" parabolic curve for a particular deal 
valuation involves establishing a limited number (e.g. three) or more 
coordinates through actual calculation, and then fitting the optimal 
curve to pass through those coordinates. 
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By way of a single-dimension simplified example/ if we consider the 
possible future values of a particular deal valuation function (d A ; 
represented graphically by the following curves over time (t) . 



Value 




(t) 

In this example if we consider this deal valuation function to be 
sensitive to a single risk factor (RFj) , say interest rates. Then, 
at a particular instant, the valuation of d*. is found by the 
operation of di(RFj) and may take the form: 

C 

Deal valuation function 



Optimal parabola 



As discussed, in order to plot each point on this deal valuation 
function would involve many thousands of test cases. Since this 
would significantly reduce the speed of operation of the method, 
rather than evaluate each point, the method makes an approximation of 
the valuation function by evaluating the deal valuation function at 
three or more values of the risk factor (RFi) and then fitting the 
optimal parabola through these points, as shown by the dotted line 
passing through coordinates A, B and C. 



Value 
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Each deal object is able to select how many coordinates to take to 
make the optimal fit parabola and is further able to choose what form 
of parabola will fit in the most suitable way. 

It should be borne in mind that the above example is highly 
simplified and that in practice the deal valuation functions and risk 
factor models will be multi- dimensional and highly complex, with each 
deal valuation function being sensitive to a plurality of such risk 
factor models. 

In this connection, a more complete analysis of the fitting of the 
optimal parabola (or quadratic polynomials) is given below. 

Fitting the Parabolic Approximation 

The key to the method is that deals of different types are able to 
select the method used to construct their optimal parabolic function. 
This parabolic function approximates the valuation function as 
closely as possible particular at the extremes of the risk factor 
where the upper-bound is most likely to be situated. 

Non-linear functions 

For non-linear products in the netted case, and all products which 
can cross the money in the non-netted case (when used without 
netting, the floor of the valuation function is used) , the following 
technique is used to fit parabolic manifolds. As already 

discussed, it is assumed that each valuation function can be 
expressed in the form: 

V(x(t),t) where x(t) is a multivariate normal distribution 
with zero mean 
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.'or risk factors that are prices, the random variable used will be 
x(t) used by the lognormal diffusion model such that: 

x(t)«$(0,<7^) 
S(t) = S(t)<e xit) 

For risk factors that are interest rates the random variable used 
will be x(t) used by the Hull-White mean reversion model such that: 



x(t)«$(0,V(« 2a< -l)/2o) 

DF{t,T) = H(tJ) • e J & r ) " 
J(t,T) = -ae' at -B(t,T) 



We seek to approximate each valuation function with a polynomial of 
the form: 



U(x, t) = a(t) + £6.(0 • x,(0 + ES^ * *. <*> • x i w 



The polynomial coefficients are generated so as to create a manifold 
that intersects the valuation function at the origin, and at the 
single dimensional upper and lower bounds of each risk factor. This 
can be regarded as a 4 discretised' Taylor expansion which tends to 
the Taylor coefficients (where calculable) as Gilt) tends to zero. 

a(t) = V(Q,t) 

C,<()= 2.(^).Q a ) 2 

Where : 

q. is a vector with the i ch element set to <T.(t) • Q a and all 
others zero 

is the normal quantile at the required level of 
confidence 
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tfhen second-order cross terms are required, they are given (for i*j) 
by: 



(n+q,- + %,t) + v(-% - q ; ,t) - v{+% - %,t) - v(-% + g,,t)) 

2.a 1 (*)-cr J (*).^ 



Where more than three points are evaluated for each risk factor., the 
process of quadratic least squares is used to fix the quadratic 
coefficients. 



Linear Functions 

For linear valuation functions in the netted case, such as forwards 
and vanilla interest rate swaps, it is possible to generate the 
coefficients of the parabolic manifold directly from the valuation 
functions themselves. Examples of direct derivation of coefficients 
are given below for a Forward Rate Agreement, and Equity Forward. 



The Forward Rate Agreement (FRA) can be written as a parabolic 
surface of two risk factors with coefficients: 



a - P 



b. = 



k = 



C 22 = 
C t2 = 



X(t) \H{t, T)-K- H(t, T + m)} 

X(t) {H^ T) - K ■ H(t,T + m)} 

X(t) ■ {H{t, T) - K ■ H(t, T + m)} / 2 

X(t) ■ (H(t, T)J(t, T)-K- H(t, T + m)J(t, T + m)) 

X{t) (H(t,T)J\t,T) - K ■ H{t,T + m)J\t,T + m))/2 

X(t) • (H(t,T)J(£T) - K • H(t,T + m)J(t,T + m)) 



K = (1 + (m/ 360) FRA) 



where x l is the risk factor for the exchange rate to base currency, 
and x 2 is the risk factor for the FRA interest rate. 



WO 03/036533 



PCT/GB01/04605 



- 14 - 



Similarly, the Equity Forward can be expressed in terms of the risk 
factors for equity price (xj deal exchange rate (x 2 ) and deal 
interest rate (x 3 ) using the following coefficients: 



o = S(t) • e- 0(r " 0 - S_FWD • H{t,T)'X(t) 

b 2 =-S_FWD-H(t,T)-X(t) 
b 3 =-S_ FWD ■ H(t, T) ■ X(t) • J{t,T) 
c u =5(i)-e-^- t »/2 
c n =-S_FWD'H{t,T).X{t)/2 
Cn=S_FWD.H(t,T).X(t).J%T)/2 
c 23 = S_FWD ■ H{t,T)-X{t) -J{t,T) 



Variations 

Some deal functions work well with a hybrid of the discretised and 
Taylor methods. For example, the value of a European equity call 
option is given by: 



V(t) = X{t)[S d (t) ■ e- D(r -° NidJ-K- DF d (/, T) • N{d 2 )] 



The value of the option in the deal's domestic currency is the 
expression within the square brackets which becomes increasingly non- 
linear as deal approaches maturity. However the exchange rate which 
converts the function back to base currency can be modelled with a 
Tayl«r expansion: 



The non-linearity of the option equation in base currency can thus be 
modelled by taking the product of the two expansions resulting in a 
closer fit polynomial. 



X(t) 



a 



c 



b 



a + bx(t) + cx 2 (t) 
X(0)-e(*'-* ;/l >' 

X{0)-e^-' l/2 > t /2 



WO 03/036533 * PCT/GB01/04605 

- 15 - 

This third operation of each deal object therefore builds parabolic 
functions which represent valuations of each deal function taking 
into account the risk factors to which that deal function is 
sensitive. 

Upper-bound of the parabolic manifold 

In order to evaluate the upper-bound of the portfolio the method uses 
the upper bound of the parabolic manifold of the portfolio as an 
approximation. The portfolio parabolic manifold is built under a 
calculation object, which in essence involves adding the parabolic 
functions of each of the deal function valuations built in the deal 
objects discussed above. 

Hence if there are 3 deals, the parabolic functions of these are 
added to arrive at a total. In practice in building the parabolic 
manifold, the coefficients of corresponding terms are added such that 
for each deal valuation function d*. 

di Vi = a x + b x . x + Ci.xr 
d 2 V 2 = a 2 + b 3 , x + c 2 ,xr 
di V 3 = a 3 + bj. x + Cj.x 2 



SUM Vj = A + B. x + C.x 2 

Where A - a 2 + a 2 + as 
B = b x + b 2 + b 3 
C » a + c 2 + cj 

This parabolic manifold will relate to a highly complex surface. 
The upper bound of this manifold is then evaluated by certain 
mathematical methods, such as aggregation, decomposition of the 
covariance matrix, transformation onto orthogonal variables, the 
application of PCA (Principal Component Analysis) on weighted 
factors, the transforming away of cross terms, and the calculation 
of distribution moments and percentiles for the approximation. These 
methods are set out below in greater detail in the context of the 
process of establishing an upper-bound for the parabolic manifold at 
a given point in time. 
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Step 1: Portfolio parabolic manifold 

As already discussed, - the portfolio is approximated by a parabolic 
manifold of normal variates with constants A(t)m Bi(t) and Cij(t) of 
the form: 

* h (0.1) 

*i«-I#<o 

k 



Step 2 : Decompose the covariance matrix 

In preparation for transforming onto an orthogonal set of risk 
factors, we find the * square root; (Cholesky factorisation) of the 
covariance matrix COV for the stochastic variables Xi. This yields a 
matrix D such that: 

COV = D-D T 

Note that the Cholesky factorisation of a covariance matrix can be 
derived from the decomposition of the correlation matrix as follows. 
Let the matrix A ~[aij] be the decomposed correlation matrix such 
that: 



AA T = 



Then the decomposition of the covariance matrix is given by [d 
a Lj ] since: 



= diag(cr)i4 • A T diag(of 

= diag(^) ^ ] diag(o) = COV 
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Step 3: Transform onto orthogonal variables 

We make the transformation onto the orthogonal vector, defined by: 

y(t) = D- l 'x(t) or x(t) = Dy{t) 

Then %(t) is multivariate normal with zero mean and unit covariance 
matrix. We may now write the portfolio approximation as: 

n(yj) = A(t)+B T -y(t)+/(t)-C-y(t) 

where : 

B = D r - B{t) 
C = D T C(t)-D 
7(0 -0(0,1) 



Step 4: Perform PCA on weighted factors 

At this stage the dimensionality of the calculation may be quite high 
- particularly if the portfolio contains a large number of equity or 
commodity derivatives. As long as these outlying assets have been 
correlated reasonably to more dominant risk factors (such as stock 
indices, the price of their trading currency etc), then movements in 
these variables is encapsulated in the dominant terms of the new 
orthogonal factors, 

To establish the principal components of the ^ factors, we take the 
absolute size of each first order term as a measure of its 
significance to the portfolio value at that timescep. 



Suppose B = 



4 



\p mJ 



then we order b. such that \b { £ b 2 



In parallel with this process, for each switch in position of the 
elements of S, we interchange both rows and columns of C. Thus if we 
switch elements 6. and bj we also interchange i th and j ch column and row 

of C. (This maintains consistency between the implicit permutation 
applied to the factor vector) . 



WO 03/036533 



- 18 - 



PCT/GB01/04605 



The elements of the orthogonal vector y are now in principal 
component order, thus in order to reduce the dimensionality of the 
calculation the first k elements may be used and the others set to 
zero. Empirically, the calculation functions satisfactorily with the 
number of principal components set between 10 and 20. 

Step 5: Transform away the cross-terms 

The elements of the C matrix still contain non-diagonal elements, so 
we make one last transformation to simplify the representation. Note 
that after PCA has been applied, this step is performed on a much 
reduced matrix, which greatly increases speed of the o(ri*) 
diagonalisation algorithm. We determine an orthogonal matrix E, and 
a diagonal matrix C(/)such that 

£(t) = E T -C{t)E and £-£ r =/ 

Then make the transformation onto the vector, z, defined by: 
L<!) = E T -M or y{t) = E-z(t) 

Note that since£-£ r =/, the vector z(C) is still (a) orthogonal and 
(b) multivariate normal with zero mean and unit covariance matrix. 
We may now write the portfolio approximation as: 

n(i('),0 = 4/) + 5(0 r 'r(0+l r W'C(/)-:(0 

whejje: 

B = E T B(t) C = £ r -C(/)£ £(0-^(0,1) 

We have thus diagonalised both the covariance matrix and the 
portfolio "Gamma* matrix C (see also Rouvinez, 1997) . 
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Step 6: Distribution moments 

We have transformed the representation of the portfolio valuation 
function into a quadratic form of independent unit normal variables. 
Leaving aside the constant term, Aft), which will be added to the 
percentile in the last stage of the calculation, we define the 
variable X as a quadratic form in normal variables: 



X = b'y + y'Cy 



A standard application of moment generating functions (see Johnson & 
Kotz, 1995) yields the following expressions for the moments (or 
cumulants) of X: 



Mean: /X = E(X) = tr(C) 

Variance : V = E((X - fxf ) = 2tr{C 2 ) + b'b 

E [{X-l*f) 8tr(C*) + 6b'Cb 



Skew: £ = 



3 3 



Kurtos 



E({X-(if) A8tr(C*) + m'C 2 b 
is: « = ; 3 = : 



More generally, the r th cumulant is given by: 

_ E((X - p.)') _ 2 {r ' i] (r - l)ltr(Cy + r\2 {r ' z) b f C {r ^b 



« r - — 

(J <T 



Although these expressions are valid for quadratic forms containing 
non-diagonal elements in the C matrix, the removal of off-diagonal 
elements greatly speeds up "the calculation of moments because (a) the 
trace {tr(C))of a diagonal matrix is just the sum of the diagonal 
elements and (b) multiplication of diagonal matrices is an o(n) 
rather than an o(n 3 ) operation. 
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:ep 7: Percentiles for the approximation 
rrom the previous step, the portfolio approximation can be written 



as: 



Where X(t) is a quadratic form of normal variables where all the 
cumulants are known. The Cornish Fisher moment correction technique 
is used to estimate the percentile (i.e. the exposure number) of this 
distribution. The q-percentile of a compound distribution of X 
involving four moments is given by: 



X a = a + a I? + i( g 2 - 1)e + ±(q 3 - 3g)« - ±(2q 3 - Sq)e* + J 



where: 



q is the percentile of the unit normal distribution. 

H, a, e, K are the mean, standard deviation, skew and kurtosis 
of X 



The q-percentile for a distribution using the first six moments is 
given by: 



+ J s (q* - 6q 2 + 3)« 5 - ^(q* - 5q 2 + 2)«yc 4 
+i(12<7*-53<r+17)« 3 

- T i 5 (2^-17g 3 + 21 g )« A 
-i(3g s -J24g 3 +29g)« 4 2 
+J r (14( 7 s -103 9 3 +107 9 )«3 2 « 4 
-^(252g 5 -1688 9 3 +1511g)«J + ... 



The exposure of the portfolio is finally given by: 



TT(t) = A(t) + X a {t) 



Note that this method functions for both netted and non-netted 
aggregations since these are determined by the first stage of 
valuation function approximation. 
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In a final operation, undertaken by the Portfolio Model object, the 
enhanced parabolic manifold, formed in the calculation object, is 
operated, i.e. it provides values for the deal valuation functions 
taking into account the controlled perturbations of the risk factors 
in the risk factor objects. 

The present invention is not to be limited in scope by the specific 
embodiment described herein. Indeed, various modifications of the 
invention will become apparent to those skilled in the art from the 
foregoing description and accompanying figure. Such modifications are 
intended to fall within the scope of the appended claims. 
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Claims 

1. A method for evaluating the credit exposure of a portfolio of 
one or more financial instruments, the method comprising :- 

establishing a deal object for the or each financial 
instrument/ the deal object comprising a representation of said 
financial instrument and a valuation function for representing how 
the value of the financial instrument is related to underlying market 
variables ; 

establishing one or more risk factor models, the or each risk 
factor model representing an underlying financial market variable 
which may affect the value of one or more of said financial 
instruments; 

establishing a deal parabolic function for representing the or 
each deal object valuation function, by operation of said the or each 
deal object valuation function on said the or each risk factor model 
to which it is sensitive; and 

summing the coefficients of each said deal parabolic function 
established at a same instant from the one or more deal objects 
represented in the portfolio in order to build a portfolio parabolic 
function which approximates the overall portfolio value for that 
instant; 

wherein establishing each said deal parabolic function involves 
evaluating a plurality of coordinates calculated from said deal 
object valuation function and then applying a parabolic curve, 
surface or multi- dimensional surface (manifold) to fit said 
coordinates, the parabolic curve or surface then representing that 
deal object valuation function. 

2. A method according to claim 1, wherein three or more 
coordinates are calculated. 

3 . A method according to claims 1 or 2 wherein the or each deal 
object valuation function is established using a deal object, each 
deal object operating to select an appropriate skeleton or template 
representation from a store of such representations of financial 
instruments, and to populate the skeleton or template with data which 
will define the nature of the financial instrument. 
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4. A method according to any preceding claim, wherein the or each 
deal object further operates to identify to which risk factor models 
the valuation function of that deal object is sensitive. 

5. A method according any preceding claim, wherein the or each 
deal object further operates to value the deal represented therein in 
relation to the risk factor models that have been identified as being 
appropriate to the valuation function of that deal object. 

6. A method according any preceding claim, wherein the deal object 
operates to apply an "optimal fit" parabolic curve to the coordinates 
evaluated from the deal object valuation function, the optimal fit 
curve being chosen on the basis that it most readily passes through 
the coordinates plotted with minimum deflection between coordinates. 

7. A method according to claim 6 wherein when more than three 
coordinates per risk factor model are evaluated, the technique of 
quadratic least squares is employed to derive the parabola 
coefficients . 

8. A method according to any preceding claim, wherein the 
coordinates to which the parabola is to be applied are chosen from 
each of the limits of the deal object valuation function (defined by 
the extremes of the relevant risk factors) and at its central region. 

9. A method according to any preceding claim, wherein templates or 
skeleton representations for each said deal object are stored in a 
suitable store. 

10. A method according any preceding claim, wherein the or each 
risk factor model is established using a risk factor object, each 
risk factor object operating to select an appropriate skeleton or 
template representation from a store of such representations of risk 
factors and to populate the skeleton or template with data which will 
define the nature of the risk factor model. 

11. A method according to any preceding claim, wherein the method 
further comprises compiling a database or store of financial and 
market data. 



WO 03/036$33 



- 24 - 



PCT/GB01/04605 



12. A method according to any preceding claim, further comprising 
the use of techniques of linear algebra and quadratic forms to 
evaluate the statistical upper-bound of the portfolio parabolic 
function , and hence, an approximation for the credit exposure of the 
portfolio. 

13. A method according to claim 12, wherein said techniques 
comprise one or more of aggregation, decomposition of the covariance 
matrix, transformation onto orthogonal variables, the application of 
PCA (Principal Component Analysis) on weighted factors, the 
transforming away of cross terms, and the calculation of distribution 
moments and percentiles for the approximation. 

14. A computer program comprising program instructions for carrying 
out the method according to any preceding claim. 

15 A computer program according to claims 14, the program being 
embodied on a record medium, computer memory, read-only memory, or 
electrical carrier signal . 

16. A method according to any one of claims 1 to 13, wherein the 
method is implemented on a computer. 

17. A system for evaluating the exposure of a portfolio of one or 
more financial instruments, the system comprising :- 

means for establishing a deal object for the or each financial 
instrument, the deal object comprising a representation of said 
financial instrument and a valuation function for representing how 
the value of the financial instrument is related to underlying market 
variables; 

means for establishing one or more risk factor models, the or 
each risk factor model representing an underlying financial market 
risk factor which may affect the value of one or more of said deal 
objects; 

means for establishing a deal parabolic function for each deal 
object, comprising means for operation of the deal object valuation 
function on said the or each risk factor model to which it is 
sensitive; and 
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means for summing each said deal parabolic function established 
at a same instant from the deal objects represented in the portfolio 
in order to build an overall portfolio value for that instant; 

wherein said means for establishing each said deal parabolic 
function comprises means for making an approximation of said each 
deal object valuation function, involving the evaluation of a 
plurality of coordinates calculated from said deal object valuation 
function and then applying a parabolic curve, surface, hyper surface 

18. A system according to claim 17, wherein the system is a 
computer system. 
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